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We consider Exclusive Double Diffractive Events (EDDE) as a powerfull tool to study 
O ! the picture of the pp interaction. Calculations of the cross-sections for the process 

^ \ p + p— >p + M + j9 are presented in the convenient form for further experimental 

1/^ ' applications. We propose measurements of t-distributions in the joint CMS-TOTEM 

experiment. It is shown that important information on the interaction region could 

be extracted from the diffractive pattern. 
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1 Itroduction 

With the first LHC run coming closer the hopes for confirmation of various theory pre- 
dictions get heated. The huge amount of works is related to the search of fundamental 
particles of the Standardt Model or its extensions (Higgs boson, Superpartners, gravitons 
and so on) and to the investigations of so called "hard" QCD processes, which corre- 
spond to very short space-time scales. " Soft" diffractive proceses take in this raw its own, 
distinctive place. 

LHC collaborations aimed at working in low and high pt regimes related to typical un- 
dulatory (diffractive) and corpuscular (point-like) behaviours of the corresponding cross- 
sections may offer a very exciting possibility to observe an interplay of both regimes [I] . 
In theory the "hard part" can be (hopefully) treated with perturbative methods whilst 
the "soft" one is definitely nonperturbative. 

Below we give several examples of such an interplay: exclusive particle production by 
diffractively scattered protons, i.e. the processes p + p — * p + M +p, where + means also a 
rapidity gap and M represents a particle or a system of particles consisting of or strongly 
coupled to the two-gluon state [2]. 

These processes are related to the dominant amplitude of exclusive and semiinclusive 
two-gluon production. Driving mechanism of the diffractive processes is the Pomeron. 
Data on the total cross-sections demands unambiguosly for the Pomeron with larger- 
than-one intercept, thereof the need to take into account the "soft" rescattering (i.e. 
"unitarisation"). 

EDDE gives us unique experimental possibilities for particle searches and investiga- 
tions of diffraction itself. This is due to several advantages of the process: a) clear 
signature of the process; b) possibility to use "missing mass method" that improve the 
mass resolution; c) background is strongly suppressed; d) spin-parity analysis of the cen- 
tral system can be done; e) interesting measurements concerning the interplay between 
"soft" and "hard" scales are possible. All these properties can be realized in common 
CMS/TOTEM detector measurements at LHC ^. 

2 Exclusive double diffraction 

The exclusive double diffractive process is related to the dominant amplitude of the ex- 
clusive two-gluon production. Driving mechanism of this processes is the Pomeron. 

To calculate an amplitude of the EDDE, we use an approach which was considered in 
detail in Ref . [2] . In the framework of this approach, the amplitude can be sketched as 
shown in Fig. 1. After the tensor contraction of the amplitudes Ti 2 with the gluon-gluon 
fusion vertex, the full "bare" amplitude Tm depicted in Fig. 1 looks like 
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Figure 1: Model for EDDE 
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with the parameters of the "hard" Pomeron trajectory, that appears to be the most 
relevant in our case, presented in Table [TJ The last factor in the r.h.s. of ([1]) is 
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where /^ = — g^ ~ q^, /i = M/2, and sq is a scale parameter of the model which is also 
used in the global fitting of the data on pp {pp) scattering for on-shell amplitudes pLj. 
The fit gives Sq — 1 GeV^. If we take into account the emission of virtual "soft" gluons, 
while prohibiting the real ones, that could fill rapidity gaps, it results in a Sudakov-like 
suppression |4[|: 
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The off-shell gluon-proton amplitudes Ti 2 are obtained in the extended unitary ap- 
proach [5]. The "hard" part of the EDDE amplitude, Fgg_,M, is the usual gluon-gluon 
fusion amplitude calculated perturbatively in the SM or in its extensions. 

The data on total cross- sect ions demand unambiguously the Pomeron with larger-than- 
one intercept, thereof the need in unitarization. The amplitude with unitary corrections. 
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Figure 2: The latest data from CDF and predictions for EDDE. 

Table 1: Phenomenological parameters of the "hard" Pomeron trajectory obtained from 
the fitting of the HERA and Tevatron data (see [Ij, [B]), and data on pp (pp) scatter- 
ing pJ.The value of the Cgp is corrected in accordance with the latest data from CDF [7], 
which is depicted in the Fig. 2 with the range of possible curves. 
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j^ , are depicted in Fig. 1. It is given by the following analytical expressions: 
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where Ai^ = Ai— g^— g^, A2T = '^2+(1t+(1t^ ^-^d the eikonal function 6pp^pp can be found 
in Ref . [T] . Left and right parts of the diagram in Fig. lb denoted by V represent "soft" re- 
scattering effects in initial and final states, i.e. multi-Pomeron exchanges. As was shown 
in [8j, these "outer" unitary corrections strongly reduce the value of the corresponding 
cross-section and change an azimuthal angle dependence. 

In the equation ([T]) we present only the Born terms from amplitudes Ti^2- It is sufficient 
for |ti^2| < 1 GeV due to fast decrease of the differential cross-section in ti^2, and the 
contribution of these corrections to the total cross-section are less than several percents. 
But when we consider the diffiactive pattern in the region of 1 < |ti,2| < 5 GeV, we have 
to take into account rescattering corrections inside the amplitudes Ti 2. In this case Is in 
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Figure 3: "Soft" survival probability 5*^ as a function of y/s for masses of the central 
system 10 GeV (solid curve), 50 GeV (dashed one) and 200 GeV (dotted one). 



the equation ([T]) changes to the following expression: 
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To calcuate differential and total cross-sections for exclusive processes we can use the 
formula 



M^ 



da 



EDDE 



\y=o 



-EDDE 



da 



Jz=0 



dW^ dy d^gg^M 

fEDDE 



rf$ 



2%6 \mV 



gg^M 
S N2(ap(0)-1) 1 



462 



LS' 



S 



2 /rf^Airf^AslT^™*'^^^ 



Jd^Aid^A2\Tj 



(14) 
(15) 
(16) 



M\ 



where da"^" ^/d^gg^M is the "hard" exclusive singlet gluon-gluon fusion cross-section and 
S^ is the so called "soft" survival probability. In this work the quantity L is called g^^g^^ 
luminocity. 

The factor 5*^ is depicted in the Fig. 3 for the systems M with quantum numbers O^"*" 
(SM Higgs boson, Radion, jet-jet). For other cases it is of the same order and can be 
calculated using the Monte-Carlo event generator EDDE [9]. 
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Figure 4: Cross-sections for the process pp — i> p + 77 + p for different kinematical cuts. 
Solid and dashed curves correspond to the pseudorapidity cuts |?7^| < 2 and |?7^| < 1. a) 
v^ = 1.8 TeV, CDF cuts for ^1,2 [T], and cut on the Ery, b) ^i = 1.8 TeV, CDF cuts for 
^1,2 [Z], and cut on the M^^; c) y/s = 14 TeV, symmetric cuts 0.0003 < ^12 < 0.1, and cut 
on the Ex-y] d) ^/s = 14 TeV, symmetric cuts 0.0003 < ^1^2 < 0.1, and cut on the M^^. 

3 Results 



First of all we would like to discuss some features of the process pp — > p + 77+p, since this 
process is the standard one to obtain the model parameters. Cross-sections for this process 
are presented in the Fig. 4. It is important to note that the cut £'77 > Ecut = Mcut/2 
is used in the major part of experimental works, that is why we have to use the same 
one in our calculations. But in some theoretical works [lOl Et > Ecut means another cut 



Myy > 2Ecut, which leads to the result, similar to the one presented in the Fig. 4b. In 
this figure cross-section for \r]^\ < 2 is about two times higher than for \ri^\ < 1. Such 
difference is only possible in the kinematics, when M^^ > 2Ecut- It follows from rather 
simple calculations. Total cross-section for the process gg —>■ 77 can be represented as [llj 
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where 'qmax is the pseudorapidity cut in the central mass frame of the diphoton system. 



C-y-y is the constant. 
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Figure 5: a) Function F{r]max); b) r]max = 2 (solid curve), r]„ 
F{ATcchj^J (dotted one), E^ut = 5 GeV. 
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is depicted in the Fig. 5a. And for the process pp — > p + 77 + p from flMl) we have 
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We are interested in the ratio of total cross-sections for different rjmax- Let us consider 
first the kinematics with cuts 



In this case 
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Since in the central mass frame of the diphoton system we have M^^ = 2ET^ch rj^, in the 
kinematics with 

^^7 = o J^ > ^c«t ' 1^7! < Vmax (22) 
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we have additional cut 
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and we should use F f min 
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max,^>-x^^^. 21;^ / instead of F{r]max) in ([17]). This new 
function is shown in the Fig. 5b. The main contribution to the integral comes from small 
masses due to fast decrease in M , but in this region we have the same cut \ri\ < ArcchT ^^ 
for different rjmax- And it is easy to get the following result 
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Even if we take as = const and L^^^^ = const, we will get the ratio 1.3, and not ~ 2 as 
in [in]. This simple example shows that we should be carefull with the kinematics during 



our calculations, since this could lead to different predictions. Now we can compare our 
result with the latest data on the exclusive 77 production from CDF 
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This prediction is higher than our calculations based on the di-jet production [7]. If we 
beleive in this result, the model parameter Cgp is at least 20% higher, than our estimations. 
The possible reason is that in 77 production we use the region of small masses for the 
normalization of our parameter (and higher masses for the di-jet production), but the 
uncertainty in the scale dependence of the cross-section is rather large (factor ~ 2). This 
result from CDF may serve a good signal for the future exclusive Higgs boson production, 
since it makes the cross-section higher by about two times. 

Now we can estimate the backgrounds for the exclusive Higgs production at LHC p]. 
Rates for these processes at the integrated luminocity 100 fb~^ are summarized in the 
table [2] (parameter Cgp = 3.7, i.e. we use more optimistic values based on the data on the 
di-jet production [7j). From this table we can obtain signal to background ratio ~ 1. More 
exact estimations will be made in the nearest future after full Monte-carlo simulations. 



Table 2: Rates for the exclusive Higgs production and different backgrounds at the inte- 



grated luminocity 100 fb and AM, 
jets with b-jets is taken to be 1%. 
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4 Diffract ive patterns 

In the general agenda of the LHC experiment diffraction often looks as an "auxiliary 
tool" for other processes such as Higgs boson and exotics searches, background supression 
and so on. Nevertheless, diffractive measurements have their own classical tasks directly 
related to the angular (or t) distributions. 

Diffractive pattern is usually characterized by the peak at small values of t, and com- 
plicated structure with dips or breaks and bumps for larger t [1]. This picture reflects the 
ondulatory properties of quantum processes as contrasted to more habitual particle-like 
behaviour, and allows us to get an information about the size and shape of the strong in- 
teraction region at large distances (i.e. directly related to confinement of the QCD colour 
fields). 
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• From the diffractive pattern we extract model independent parameters of the inter- 
action region such as the t-slope which is -R^/2, with R the transverse radius of the 
interaction region. 

• We can also estimate the longitudinal size of the interaction region [13] : 

Axl > ^ (26) 

The longitudinal interaction range is somehow "hidden" in the amplitude but it is 
this range that is responsible for the "absorption strength". A rough analogue is 
the known expression for the radiation absorption in media which critically depends 
on the thickness of the absorber. 

• The very presence of dips is the signal of the quantum interference of hadronic 
waves. 

• The depth of dips is determined by the real part of the scattering amplitude 

According to the data from SppS and Tevatron the transverse radius of the interaction 
region is of order of 1.2 fm ~ 1.5 < rem >• The longitudinal size can be estimated from 
the second inequality (126|) and is of order of 2800 fm. 

Diffractive pattern moves due to changes in kinematical parameters like the energy 
of the interaction or an additional hard scale. This motion reflects the dynamics of the 
process. The increase of the t-slope with energy reflects the growth of the interaction 
radius. At fixed collision energy the diffractive pattern is fixed as well. 

However if we have in our disposal an additional hard scale we can operate the diffrac- 
tive pattern adjusting this hard scale at our will and making, e.g., the interaction region 
larger or smaller. 

Hard scale is related to small distances and, from the simple optical point of view, 
the pattern should move towards large values of -t with the increase of the hard scale. 
HERA provides an excellent opportunity to observe the influence of a hard scale (Q^) on 
the diffractive pattern: the slope decreases with Q^ in exclusive vector meson or photon 
production or for mesons containing heavy quarks [J/"^) as contrasted to those composed 
of the light quarks (see Fig. 6) [H]. 

The decrease of the slope with Q^ in electroproduction was predicted qualitatively 
in Ref. [15]: J.D. Bjorken argued that the decrease of the slope would be bounded from 
below by the size of nucleon [16] . The latter feature seems to be violated in the HERA 
data [17]. We have to mention that the presence of a high-mass particle in the final 
state does not always lead to the phenomena described above. For example, hadronic 
resonances with large masses have large size due to intrinsic motion of constituents, and 
can not be considered as a hard probe. In this case we have inverse dynamics of the 
pattern [18]. This certainly is not the case for the processes considered below as they are 
related exclusively to short-distance probes, i.e. "high mass" means always "high Et". 

The diffractive pattern for the process p + p-^p + jj+p a.s predicted on the basis of 
ref. [2] is displayed on the Fig. 7 where ^^ means the exclusive differential cross section 
with all final variables integrated except one of the proton transverse momenta (— t) and 
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Figure 6: The slope b, as a function of Q^ + My, compared to other ZEUS and HI results. 




12 3 4 5 



-t,GeV' 



Figure 7: Normalised cross section for exclusive dijet production as a function of t for 
Mx = 30 GeV (the solid and long-dashed curves correspond to the LHC and TEVATRON 
energies, respectively) and Mx = 200 GeV (the dotted and short-dashed curves corre- 
spond to the LHC and TEVATRON energies, respectively). The left curve corresponds to 
the elastic scattering at the LHC. 

the value of the central mass (M ~ 2Et). With two exclusive high-Ex jets the expected 
dips will reflect the elastic scattering of the protons off the hard gluon. Their positions are 
shifted to the right in comparison with the proton-proton elastic scattering, as depicted 
on the Fig. 7. Such a shift is a clear signal of the short-distance scale due to jets. 

Measurements of t-distributions and their dynamics in the exclusive central diffraction 
could be used for the proposed investigations. To obtain the detailed diffractive pattern 
with dips for 1 GeV^ < — t < 5 GeV^ we need at least 10^ events for fixed (or falling 
within the small enough range of values) masses of the central system and t-resolution 
less than 10% in this region. At high luminosities the use of the missing mass method 
is limited below by central masses above 30 GeV because of the acceptance limitations 
and the absence of resonances with high rates in this region. That is why the only way 
is to use exclusive or Semi- Inclusive (exclusive+" soft" radiation in the central rapidity 
region) dijet production. The best case is the measurements at the nominal luminosity 
at P* = 0.5. Results are summarized in the table [31 
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Table 3: Rates for exclusive and semi-inclusive {\risoft\ < 5) double diffractive dijet pro- 
duction for luminosity 10'^^ cm~^ s 
central system, Mx- 



33 ^v>-, 2 ^ 1 £qj, different intervals of the invariant mass of the 
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